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We calculate one loop corrections to the domain- wall quark propagator in QCD. We show how the wave function 
is renormalized in this theory. Especially we are interested in the behavior of the massless fermion mode, which 
exists near the domain wall at the tree level. We show that this massless mode is stable against the quantum 
correction. 



1. Introduction 

The theory of domain wall fermion was origi- 
nally introduced to treat the chiral gauge theory 
. Apart from this expectation the domain wall 
fermion is regarded as a suitable formulation to 
treat the vector like massless QCD This is be- 
cause of the following great advantages compar- 
ing with the ordinary Wilson or Kogut-Susskind 
fermion; (i)The number of flavors is not fixed. 
(ii)The renormalization of mass is multiplicative. 
In other words, if a massless mode exists in the 
tree level it is stable against the quantum correc- 
tion. The first one is almost trivial but the latter 
is not. It has been understood from an intuitive 
discussion and a numerical simulation [^,^. The 
aim of this paper is to confirm the latter nature, 
especially the stability of the massless mode by 
lattice perturbation theory. 

In this paper we set the lattice spacing a = 
1 and take the SU{N) gauge group with second 
Casimir C2. We adopt the domain wall fermion 
of Shamir type ||] to describe massless quarks. 
This type of domain wall fermion is a variation 
of the Wilson fermion with sufficiently many Ns 
flavors and special form of mass matrix. The only 
difference from the Wilson fermion action is the 
fermion bilinear term. In the perturbation theory 
the gluon propagator and the gauge interaction 
terms are exactly same as those in the ordinary 
Wilson fermion perturbation theory with Ns 
flavors. The peculiar thing to the domain wall 
fermion is the fermion propagator which is given 



as an inverse of the fermion bilinear term. 
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where P± is a projection operator; P± = (1 ± 
75)72. Our mass matrix W^{p) is given as 
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In spite of the Dirac mass M in the mass matrix, 
the above action describes one massless fermion 
and Ns — 1 excited modes at the momentum 
region by virtue of this matrix form. 

This means that M is not a physical quark mass 
but it is rather a cut-off order unphysical mass 
like Wilson mass. M plays an important role 
as a parameter of the theory in order to satisfy 
iW^lPp 1' which is needed for the exis- 

tence of a massless mode. 

In the following we will see this massless mode 
is stable against the quantum correction. 

2. Massless Mode 

Before calculating the one loop correction we 
show how the massless mode is guaranteed in the 
domain wall fermion system. 
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The easiest way to see the massless mode is to 
diagonahze the fermion bihnear term. The bihn- 
ear term can always be diagonaUzed by rotating 
the right and left mode of the Dirac fermion in- 
dependently by two unitary matrices U and V, 



(5) 



These two unitary matrices U and V should be 
chosen to diagonalize two different hermite mass 
squared matrices W^W~ and W~W^. For ex- 
ample the mass squared eigenvalues are following 
at M = 0.8 and Ns = 20, 



(1.43516, 



0.645049, 10- 



-18\ 



(6) 



The important point is that there is only one 
massless mode and other modes are excited ones 
with cut-off order mass. By virtue of this nature 
the domain wall fermion can be regarded as one 
massless Dirac fermion system with extra excited 
modes decoupled. 

This property of the mass eigenvalues can be 
understood from the discussion of the eigenstate 
equation of the mass squared matrix W~W^, 



{W-W+)st^l^M^C 
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where and Mf are the i-th eigenstate and eigen- 
value. By substituting mass matrix and (|^). 
we have three types of equations depending on s. 



- W {C+i + C-i) + {l + W'- Mf) C = (8) 
-WQ + {W^-M^)(l^O (9) 
-WCns-i + {l + W'- Mf) Cns = (10) 

From the equation ^ which is valid for 2 < s < 
Ns — 1 we have two kinds of solutions depending 
on the eigenvalue. When the eigenvalue is small 
Mf < (1 — WY we have exponential solution, 



C(s)=Ae±"", cosh a = 



l + W^ - Mf 
2W 
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When the eigenvalue is in the region (1 — WY < 
Mf < (1 — W)"^ we have an oscillating solution 



C(s) = Ae 



±ias 
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And when the eigenvalue is larger; Adf > (1 — VF)^ 
we again have an exponential solution. 



Besides the general equation (|8|) the solutions 
have to satisfy the boundary conditions given by 
the equations (^ and (|ic|), 
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Here we assume that Ns is infinitely large. First, 
a single dumping solution can satisfy the first 
boundary condition with a dumping factor W. 
This is nothing but a zero mode solution with 
Mf — 0. Other dumping solutions cannot sat- 
isfy this condition. On the other hand oscillating 
ones can always satisfy ( p^ ) with suitable Mf . In- 
deed we can solve the conditions in this case. As 
a result we have one dumping solution with zero 
eigenvalue and oscillating solutions with eigenval- 
ues in the region {l-W^ < Mf < (l + W)^. For 
example the tree level result at M — 0.8 has one 
small eigenvalue and Ns — 1 eigenvalues in the re- 
gion; (1 - W)^ = 0.64 < Mf <{! + Wf = 1.44. 
This is a good coincidence with the above discus- 
sion. 

3. One Loop Calculation 

Now we calculate the one loop correction to the 
fermion propagator, which is given by a tadpole 
diagram with a gluon loop and a half-circle dia- 
gram, in which both gluon and fermion runs. We 
are interested in stability of the massless mode 
given in the — > limit. We only require the 
leading terms in to see the massless mode. 

The contribution from the tadpole diagram is 



^tadp 



olc 



2 (5. 
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with numerical factor T = 0.15461. The first 
term in (|l^) is finite, and the second term is cut- 
off order. We can see that 5]*'"^?°''^ is diagonal in 
fiavor space. 

The correction from the half circle diagram 
cannot be calculated analytically because of its 
complicated dependence on the flavor index s, 
which is brought by the fermion propagator. We 
will evaluate the form of S^f ^ according to by 
separating the loop momentum into two region. 
The logarithmically divergent part can be cal- 
culated analytically and other linearly divergent 
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and finite terms are given by integrating the loop 
momentum numerically, where we let Ns — 20 
for the calculation. And finally the correction is 
written in a simple form, 



s'^'^'f = - h 



log 



finite 



(16) 



where finite term /g^j^j, and linearly divergent 
are complicated functions of s,t and M. 

The logarithmic divergence appears only in the 
wave function and is localized near the boundary 

^4 = I^5'^2Af(2 - Af)(l - Af)^+*-2 

X (^ln(7r2) + i - lnp2^ . (17) 

In the diagonal basis of (||) it is exactly localized 
in the boundary. This means that the logarith- 
mic divergence can be rcnormalizcd into the zero 
mode wave function. On the other hand, the lin- 
ear divergence in the mass term should be treated 
as an additive quantum correction to the mass 
matrix. This is because the tree level mass ma- 
trix is cut-off order and a correction of the same 
order is always allowed. 

In the following we will see whether the mass- 
less mode is preserved against the correction to 
the mass matrix. The contribution from the tad- 
pole diagram is to modify the mass parameter to 
M - 2g'^C2T. Although M± has nontrivial flavor 
dependence, it is much smaller (~ 18% at most) 
than that of the tadpole diagram and it turns out 
to be a similar form to the tree level mass matrix; 
minus values along the diagonal line and plus ones 
along next to the diagonal. Differently from the 
tree level mass matrix (||), it's off-diagonal ele- 
ments are non-zero. But they are very small com- 
paring the tadpole contribution (~ 7% at most) 
and dump exponentially along the off-diagonal 
line with dumping factor 0.12 aX M — 0.8. 

By virtue of the similarity to the tree level mass 
matrix we may expect that the good nature of the 
domain wall fermion is preserved at the one loop 
level. Diagonalizing the effective mass matrix nu- 
merically we have following sequence of the mass 
squared eigenvalues for M — 0.8, 



(2.36322, 



0.306761, 10""). 



(18) 



Here again we have a single zero eigenvalue and 
Ns — 1 cut-off order ones. Corresponding to them 
we also have one dumping solution and Ns — 1 
oscillating ones. 

From this peculiar distribution of eigenvalues 
and eigenstates we may expect that the same 
mechanism works as at the tree level and the 
above nature is guaranteed. The one loop ef- 
fective mass matrix is well approximated by the 
same form of the matrix (^). Actually with the 
approximated value of W = 0.5141 the numeri- 
cally solved exact results are well reproduced as 
an eigenvalue problem of this simple matrix. The 
above W is about 7% smaller than the dump- 
ing factor of the exact zero mode solution, which 
should be coincide in our approximation. This 
difference may be the effect of the off-diagonal 
parts which we neglected. 

Our conclusion is that the massless mode is 
preserved at one loop level because the mass ma- 
trix has approximately the same form as at the 
tree level. This fact also guarantees the other 
modes have cut-off order mass. This approxi- 
mation is supported by the following reasons; (i) 
main quantum correction comes from the tadpole 
diagram, which is diagonal, (ii) the off-diagonal 
parts are small and dump exponentially off the 
diagonal line, (iii) numerically solved mass eigen- 
values and eigenstates are well reproduced by this 
approximated form of the mass matrix. For our 
approximation to work it is essentially important 
that the off-diagonal parts dump exponentially. 

At last we comment that the finite parts in the 
wave function subtracted with the logarithmic di- 
vergence can be diagonalized simultaneously. At 
M = 0.8 the Z-factor of the zero mode is given as 

Z ^1- 0.021660252C2. 
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